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We report the results from a computational study of the self-assembly of amphiphilic ditethered 
nanospheres using molecular simulation. We explore the phase behavior as a function of nanosphere 
diameter, interaction strength, and directionality of the tether-tether interactions. We predict the 
formation of seven distinct ordered phases. We compare these structures with those observed in 
linear and star triblock copolymer systems. 
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I. INTRODUCTION 

The ability to create ordered structures on the 
nanoscale has important implications for the fabrication 
of novel electronic [H [2] and photonic [3j 0] materials. 
To this extent, the literature contains a variety of ap- 
proaches to the assembly process. For example, dipolc 
moments have been exploited to induce nanoparticles to 
form sheets [5j [6] and wires [5j |7j [8], binary nanopar- 
ticle superlattices reminiscent of atomic crystals have 
by created by controlling parameters such as particle 
size ratio and charge ratio [H E] , and block copolymers 
have been used as templates to create high conductiv- 
ity wires [5] and other ordered arrays of nanoparticles 
[TU] . The literature also contains many examples of the 
use of polymer-tethered nanoparticles as a means to self- 
assemble ordered structures pU OH H3 LH HSJ ESI El 

Tethered nanoparticles are hybrid nanoparticle-polymer 
building blocks where nanoparticles are bonded to im- 
miscible polymer tethers creating amphiphilic building 
blocks The immiscibility between the nanoparti- 

cle and polymer tether facilitates microphase separation 
into bulk periodic structures similar to those observed in 
block copolymers, but with additional ordering arising 
from the nanoparticle shape [HllI1Zlirai2Tll2ini3tj] . 

In reference [17] we studied systems of ditethered 
nanospheres (DTNS) where two chemically distinct 
polymer tethers were attached to the surface of a 
nanosphere, creating a building block with multiple levels 
of anisotropy [32] • The DTNS building block is of partic- 
ular interest as the synthesis of nanospheres functional- 
ized with diametrically opposed tethers has recently been 
demonstrated, resulting in the formation of wires under 
dilute conditions |21[ 122] and preliminary indications of 
cylinders [33], as predicted in earlier simulations 1 1 J . 
Just as mono-tethered nanoparticles have been shown to 
form phases similar to diblock copolymers [TT] [T2 [TS] 
and single-tailed surfactants [HI [14], DTNS have simi- 
lar analogy with double-tailed surfactants and triblock 
copolymers [17]. DTNS can be loosely thought of as a 
nanoparticle equivalent of an ABC triblock copolymer 
where the center block of a triblock has been replaced 
by a nanosphere. In our previous study of DTNS, we 



found a variety of structures analogous to triblock copoly- 
mers, including the alternating double diamond network 
and alternating tetragonal cylinders [17]. We also found 
two novel structure not seen in pure triblock copoly- 
mer melts, namely NaCl ordered spherical micelles with 
a complementary simple cubic network of nanospheres 
and zincblende ordered spherical micelles with a comple- 
mentary diamond network of nanospheres |17j . In this 
work we explore the phase behavior of the DTNS sys- 
tem as a function of nanosphere diameter, immiscibility, 
and directionality of the tether-tether interactions to bet- 
ter understand the structural behavior, expanding on the 
previous study in reference |17j . 

In section[TT]we discuss the simulation model, method, 
and analysis methods used in this work. In section [ITl] we 
present the phase behavior as a function of nanoparticle 
diameter, immiscibility, and directionality of the tether 
interactions. In section HVl we discuss these results and 
in section |V] provide concluding remarks. 




FIG. 1: Schematic of the DTNS building block. The blue 
colored tether is labeled as A, the nanoparticle as B, and the 
red tether as C. 



II. MODEL AND METHODS 

A. Simulation Model and Method 

In the DTNS system, self-assembly is driven by immis- 
cibility and mitigated by (1) the angle between the teth- 
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ers, (2) the size and shape of the nanoparticle, and (3) 
the interaction strength between like species. We utilize a 
minimal model that captures these important effects and 
use empirical pair potentials that have been successful in 
the study of the phase behavior of block copolymers, sur- 
factants, and colloids, and in previous studies of tethered 
nanoparticles. 

We model polymer tethers as bead spring chains, where 
individual beads of diameter a are connected via finitely 
extensible non-linear elastic (FENE) springs (34] , with 
the maximum allowable separation set to Rq — 1.5a and 
the spring constant set to k = 30 [33]. Nanospheres are 
modeled as beads of diameter D connected to the poly- 
mer tethers also via FENE springs. The planar angle, 
9, between tethers at the surface of the nanoparticle is 
controlled by the use of a harmonic spring, with k = 30 
and i?o = (D + a)sin(8/2). A schematic of the building 
block is shown in figure [l] We fix the length of the two 
tethers at 4 beads each, and vary nanosphere diameter, 
D = 1.5<t, 2.0a, or 2.5a and planar angle, 9 = 30° to 
180°. 

We use the Brownian dynamics (BD) simulation 
method [33] to realize long time scales and large systems 
required to self-assemble complex mesophases. We per- 
form simulations under melt-like conditions where like 
species attract and unlike species do not attract. To 
model the attraction between like species, we use the 
Lennard- Jones potential (LJ), which induces demixing 
below a certain critical temperature. The LJ potential is 
given by, 

Ul ] = f 4e ( (r-a) 12 ~ W=w) + Ushift ' r ~ a < Tcut 
1.0, r-a>r cut 

(1) 

where e is the attractive well depth, U shift is the en- 
ergy at the cutoff, a is the parameter used to shift the 
interaction to the surface of the nanoparticles to prop- 
erly account for excluded volume, and r cut is the cutoff 
of the potential with respect to a. a is determined by 
calculating a = (Di — a)/2 + (Dj — a)/2, where Di and 
Dj are the diameters of the two interacting particles. 
Note that for interactions between tethers of like species, 
where Di = Dj = a, a = 0. The cutoff of the LJ poten- 
tial is set to r cut = 2.5a. Species of different type interact 
via a purely repulsive Weeks-Chandler- Andersen (WCA) 
soft-sphere potential to account for excluded volume in- 
teractions. The WCA potential can be described by the 
LJ equation (cqn. [I]), with U s h%ft = e and r cut — 2 1 / 6 <7 

Volume fraction, ip, is defined as the ratio of excluded 
volume of the beads to the system volume. The degree 
of immiscibility is determined by the reciprocal temper- 
ature, 1/T* = e/khT. In this study, the potentials are 
chosen such that they capture the overall nanoparticlc- 
tether and tether-tether immiscibility, the geometry of 
the nanoparticle, and the angle between the tethers. 
Changes to the phase behavior are likely if the individ- 
ual interaction strenghts are changed asymmetrically, but 



such changes are not within the scope of this study. 

We perform simulations at </> = 0.45 utilizing the fol- 
lowing simulation procedure. For a given 9 and D we 
start from a high-temperature, disordered equilibrated 
state and incrementally cool the system, allowing the po- 
tential energy to equilibrate for several million time steps 
at each T* before cooling again. For each 9 we perform 
at least two different cooling sequences in all cases to as- 
certain the path independence of the observed phase and 
perform simulations at different system sizes to eliminate 
artifacts due to box size effects. 



B. Method for calculating micelle centers-of-mass 

To analyze structures composed of spherical micelles, 
it is often necessary to locate the center-of-mass of each 
micelle. To approximate the center-of-mass we modify 
an image processing technique developed by Crocker and 
Grier that is typically used to identify the center of col- 
loidal particles from microscopy data [36] [37]. In this 
method, we start by constructing a density profile of the 
system by creating a 3-d grid of cubic cells and calcu- 
lating the number of particles in each cell. To increase 
our accuracy, we use approximately 50 snapshots spaced 
50 timesteps apart to average over the time-dependent 
shape of the micelles, but within a small enough time 
window that the bulk structure is still correlated. We 
then multiply the density field by a 3-D Gaussian func- 
tion (i.e. apply a 3-D Gaussian filter) [3S] [37] |3S] where 
the diameter of the Gaussian corresponds to the approx- 
imate diameter of the micelles. This weights the center 
of a spherical object greater than the edges. We then 
calculate the cell with the highest weighted value in each 
region, where a region is defined by the approximate di- 
ameter of the micelles. This calculation produces a set of 
cells that correspond to the centers of mass of the micelles 
in the system. To avoid artifacts associated with the res- 
olution of the grid, we vary our cell size between 0.25 and 
0.5 . We additionally calculate the centroid around each 
center-of-mass cell and then average the values to arrive 
at the centers used in our analysis. 



C. Mesostructure identification 

To identify the mesophases formed by the DTNS we 
utilize a combination of visual inspection, calculation of 
the structure factor, S(q) [351 [3D], and construction of 
the bond order diagram (BOD). The structure factor is 
a mathematical description of how a structure scatters 
incident radiation due the arrangement of material. The 
structure factor produces a set of strong peaks whose ra- 
tio can be used to identify specific geometric structures 
[40] . The BOD is constructed by taking the directions of 
all vectors drawn from a particle or micelle to the nearest 
neighboring particles/micelles. These vectors are then 
projected on the surface of a sphere, creating an "aver- 
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age" picture of the orientational order in the system. Sys- 
tems that have highly correlated nearest neighbor direc- 
tions (e.g. bulk crystalline materials) will show distinct 
groups of points on the surface of the sphere. Disordered 
systems, such as liquids, will appear as points randomly 
distributed on the surface of the sphere with no clearly 
correlated points. For the calculation of these quantities 
we average over several independent simulation runs and 
up to 100 time slices. 



III. PHASE BEHAVIOR AS A FUNCTION OF 9, 
D AND IMMISCIBILITY 



In this section we explore the phase behavior of diteth- 
ered nanospheres under melt-like conditions as a function 
of planar angle, immiscibility, and nanoparticlc diameter. 
Note that the length of each tether is fixed at 4 beads. 
The overall phase behavior, grouped by nanosphere di- 
ameter D, is summarized in the "phase diagrams" shown 
in figures [2^,-c; over 500 individual statepoints are in- 
cluded in these diagrams. The individual structures 
found in figures |2]i-c are discussed in detail in the sub- 
sections that follow. Throughout this chapter, figures are 
color coded as shown in figure [T] A tethers and the ag- 
gregates they form are colored blue, B nanoparticles and 
the aggregates they form are white/gray, and C tethers 
and the aggregates they form are red. 



A. Lamellar sheets/nanosphere monolayer (L/M) 

A simulation snapshot of a lamellar structure (L/M) 
for L>=1.5er, 6>=150° and 1/T*= 0.8 is shown in figure [3] 
Nanospheres form monolayers sandwiched between alter- 
nating layers of tethers. The spacing between the layers 
varies slightly with 9 with the largest spacing occurring 
for systems where 0=180°. For lamellae, the spacing be- 
tween layers (i.e. periodicity) can be calculated using 
L spacing = (211/(7*), where q* is the modulus of the wave 
vector at which the first principal maximum in S(q) is 
located gHHS]. We find that for D=1.5a and 9 = 90°, 

3.9o\ 



spacing 



150°, and 180°, the average spacing is L 
4.2er, and 4.3(7, respectively; the spacing does not appear 
to depend on 1/T* given we are in the ordered regime. 



B. Alternating tetragonally ordered 
cylinders/nanosphere cylinders (TC/C) or 
tetragonal mesh (TC/T) 

Each nanosphere diameter investigated demonstrates 
a region of cylindrical structures. In all cases, we 
find a cross-sectional checkerboard pattern of alternat- 
ing tetragonally ordered cylinders formed by the tethers. 
However, the arrangement of the nanospheres is different 
when we consider different nanosphere diameters. For 
D— 1.5a nanospheres organize into cylinders (TC/C), as 
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FIG. 2: Phase diagrams as a function of planar angle, 8, and 
1/T*, grouped by D . The structure labels are as follows: al- 
ternating tetragonally ordered cylinders formed by the tethers 
where the nanospheres also form cylinders (TC/C); alternat- 
ing tetragonally ordered cylinders formed by the tethers where 
the nanospheres form a tetragonal mesh (TC/T); alternating 
tetragonally ordered cylinders formed by the tethers where 
the nanospheres form cylinders or a tetragonal mesh (TC); 
alternating double gyroid structure within a nanosphere ma- 
trix (AG); alternating double diamond structure within a 
nanosphere matrix (AD); lamellar sheets where nanospheres 
form a monolayer (L/M); NaCl ordered spherical micelles 
where nanospheres form a complementary simple cubic net- 
work (NaCI/SC); and ZnS (i.e. binary diamond) ordered 
spherical micelles where nanospheres form a complementary 
diamond network (ZnS/D). Phase boundaries are approxi- 
mate and drawn to guide the eye. c is adapted from reference 



FIG. 3: Lamella structure formed for D=1.5<r, 0—150° and 
1/T*= 0.8 showing a monolayer of nanospheres. The spacing 
between the nanosphere layers is ~4.2<r. 



shown in figure [3Ji. In contrast, for D=2.ha nanospheres 
organize into a tetragonal mesh that separates the al- 
ternating cylinders (TC/T), as shown in figure Bp. In 
both cases, the structures can be described by the[8,8,4] 
Archimedean tiling constructed of octagons and squares, 
shown in figure [4J; for D—1.5a and figure|4]i for D=2.5a. 
We should note that for D=1.5a tethers occupy the oc- 
tagonal tiles and nanospheres occupy the square tiles, 
whereas for D=2.5a both tethers and nanospheres oc- 
cupy the octagonal tiles and only nanospheres occupy 
the square tiles. 

For D=2.0a we see a both tetragonal and cylindri- 
cal morphologies of nanospheres. As the value of 1/T* 
increases (i.e. we cool the system), there is a gradual 
transition from a tetragonal mesh of nanospheres to dis- 
tinct cylinders formed by the nanospheres. Conceptually, 
the tetragonal phase can be thought of as the cylinder 
phase with connections between the cylinders. As the 
strength of the attraction between the nanospheres in- 
creases (i.e. we cool the system), these connections are 
minimized forming more compact, cylindrical nanosphere 
domains. 



C. Alternating gyroid/nanosphere matrix (AG) 

A simulation snapshot of the alternating gyroid struc- 
ture (AG) is shown for D=2.0a, 0=180°, and 1/T* = 0.4 
in figure |5^. The AG structure consists of two chemically 
distinct gyroid networks, one network formed by the A 
tethers and one formed by the C tethers, separated by a 
matrix of nanospheres. These networks are composed of 
a series of interconnected tubes, where three tubes meet 
at each node. The AG structure has the same Ia3d space 
group of the double gyroid (DG) [16j |43] , however each 
gyroid network is chemically distinct in the AG phase. 
This structure was identified both visually and by the 
use of the structure factor, S(q), finding characteristic 




FIG. 4: (a) Simulation snapshot of the end view of the TC/C 
phase for D=1.5a, 9 = 60°, and 1/T* = 0.8. (b) Simulation 
snapshot of the end view of the TC/T phase for D=2.5a, 9 = 
90°, and 1/T* = 0.8. (c) [8,8,4] Archimedean tiling overlayed 
on the simulation snapshot in (a), (d) [8,8,4] Archimedean 
tiling overlayed on the simulation snapshot in (b). 



peaks in ratio y/3 : y4 : y/l0 : vlT as expected [H] and 
shown in figure [5p; for ease of viewing the x-axis was 
scaled by 1 .47 such that the numerical values correspond 
to the values in the characteristic ratio (i.e. the first peak 
occurs at a numerical value of y/S). 



D. Alternating diamond/nanosphere matrix (AD) 

The alternating diamond phase (AD), shown in fig- 
ure [5ji for D=2.5a, 9 =180° and 1/T*=0.8, consists of 
two chemically distinct, interpenetrating diamond net- 
works, one formed by the A tethers and one formed by 
the C tethers, separated by a matrix of nanospheres. 
This phase could also be classified as the double diamond 
structure if we ignore the chemical specificity of the two 
diamond networks |44j . Each diamond network is com- 
posed of cylindrical tubes, where four tubes connect at 
a node in a tetrahcdral arrangement. This structure was 
identified visually and by calculating S(q). S(q) for the 
alternating diamond phase is plotted in figure [5jl, show- 
ing characteristic peaks with ratio v2 : y/3 : v4 : as 
expected [H]; for ease of viewing the x-axis was scaled 
by 1.34, such that the numerical values correspond to 
the values in the characteristic ratio (i.e. the first peak 
occurs at a numerical value of y/2). 
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FIG. 5: (a) Isosurfaces of the two tether domains, which self- 
assembled into the AG structure showing 8 unit cells with a 
unit cell size of 13a; nanospheres are rendered at half their 
true size, (b) Structure factor of the AG phase. Data in 
(a-b) correspond to D=2.0cr, 6 = 180° and 1/T* = 0.4. (c) 
Isosurfaces of the two tether domains, which self-assembled 
into the AD structure showing 27 unit cells with a unit cell size 
of 10a; nanospheres are rendered at half their true size, (d) 
Structure factor of the AD phase. Data in (c-d) correspond 
to D=2.5cr, 6 = 180° and 1/T* = 0.8 and are adapted from 
reference |17j . 



E. NaCl ordered spherical micelles/nanosphere 
simple cubic network (NaCl/SC) 

For D— 2.0a and 2.5a the tethers form spherical mi- 
celles that order into an NaCl lattice with a comple- 
mentary simple cubic network formed by nanoparticles 
(NaCl/SC). In figure |6^i we plot a simulation snapshot of 
the centers of mass of the micelles formed by the tethers 
that order into a NaCl lattice for D=2.5a, 8—60°, and 
1/T*=0.8. We can see that the structure clearly demon- 
strates alternating chemical specificity. Figure [BJd shows 
the bond order diagram (BOD) for the centers of mass of 
the micelles, where we ignore chemical specificity of the 
micelles and simply calculate the BOD for nearest neigh- 
bor micelles. The resulting BOD shows a simple cubic 
arrangement of the micelles, which corresponds to the 
BOD of a perfect NaCl system, shown as lines in figure 

The nanoparticles fill the space between the micelles 
arranging into a simple cubic network; an isosurface of 
the nanoparticles is shown in figure [6]i. The nodes of the 
simple cubic structure each have six connections points; 
a single node is shown in figure |6j;. The highest den- 
sity of nanoparticles (i.e. the node of the SC network) 
resides in the center of an eight-particle NaCl unit cell, 
as shown in figure [6J3 - this high density location corre- 
sponds to the placement of the central particle in a BCC 



FIG. 6: (a) Centers of mass of the NaCl ordered micelles 
formed by the tethers; the unit cell size is approximately 6.5a. 
A perfect NaCl unit cell is inset, (b) BOD of the micelles' 
center of mass; the BOD of a perfect NaCl structure is shown 
as lines, (c) Isosurface of the nanoparticles showing a sim- 
ple cubic network arrangement, (d) Eight particle unit cell 
of NaCl formed by the micelles in the system, with the node 
of the nanoparticles network at the interstitial shown as an 
isosurface. (e) 8 particle unit cell of NaCl formed by the 
micelles in the system with the nodes of the nanoparticle net- 
work drawn as white spheres showing the AlCu2Mn structure. 
All data for D=2.5a, 6 = 60° and 1/T*=0.8. Adapted from 
reference |17j . 



lattice. We approximate the location of the nodes using 
the micelle center of mass calculation; figure shows an 
eight-particle unit cell of NaCl extracted from our system 
with the locations of the nanoparticle nodes rendered as 
spheres. The overall phase corresponds to the AlCu2Mn 
structure, also known as the Heusler (L 2 1) phase [35] . 
The AlCu2Mn structure is a three-component analog of 
BCC (note, CsCl is a two-component analog to BCC). 



F. Zincblende ordered spherical 
micelles/nanosphere diamond network (ZnS/D) 

For D=2.5a we find a binary mixture of spherical mi- 
celles formed by the tethers that order into a zincblende 
(binary diamond) lattice with a complementary diamond 
network formed by the nanospheres (ZnS/D). In figure 
[7^1, we show a simulation snapshot of the centers of mass 
of the micelles formed by the tethers in the zincblende 
structure. This structure is the two-component analog of 
the diamond lattice [IS]. Figures |7|3-d show the BODs of 
the micelles' centers of mass for the zincblende structure. 
We split the BOD into two separate diagrams, since the 
diamond phase possesses two bond configurations that 
are 60° rotations of each other. We plot clusters where a 
"blue" particle (i.e. micelle formed by the A portion) is 
at the center surrounded by "red" particles (i.e. micelles 
formed by the C portion), and a second diagram where 
clusters have a "red" particle at the center surrounded by 
"blue" particles; these cluster definitions properly group 
the data by orientation of the tetrahedrons. Both BODs 
in figure [TJd and c show clear tetrahedral arrangements; 
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we plot the BODs for an ideal zincblende structure as 
lines in both plots, showing good agreement with our 
simulation results. The blue centered and red centered 
tetrahedrons have complementary orientations (i.e. rota- 
tions of 60°), as shown in figure |7ji. 
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FIG. 7: (a) Centers of mass of the micelles formed by tethers 
that order into the ZnS lattice showing 27 unit cells with a 
unit cell size of approximately 8<r. A perfect diamond lattice 
is shown in the inset, (b-c) BOD for nearest neighbors in the 
ZnS structure for the two different tetrahedral arrangements, 
(d) The combination of figures (b) and (c). (e) Diamond 
network formed by the nanoparticles in the Zns/D phase, (f) 
Diamond network formed by the nanoparticles in the Zns/D 
phase, (f) Average S(q) for the nanoparticle network, show- 
ing strong peaks with ratios v3 : VS : VTT : \/l6. Note that 
q was scaled such that the first peak corresponds to v3 for 
ease of viewing. All DTNS data at D=2.5ct, 9 = 30° and 
1/T* = 0.8. Adapted from reference [17] . 



Within this phase, the nanoparticles organize into a 
diamond network that is woven into the micellar lat- 
tice. An isosurface of the nanoparticle diamond net- 
work structure is shown in figure [7£. In figure [Tf we 
plot S(q) for the nanoparticles, finding strong peaks at 
ratios \/3 : a/8 : \/Tf : \/l6 as expected for diamond 
[10]; note that q was scaled such that the first peak cor- 
responds to a/3 for ease of viewing. The overall ZnS/D 
phase is composed of two interwoven diamond structures 
where one network is formed by nanoparticles and the 
other consists of a binary lattice of spherical micelles. 



IV. DISCUSSION 

The overall phase behavior as a function of nanosphere 
diameter can be better understood by constructing a plot 
of D verses 9 for a fixed 1/T*. In figure [8] we present the 
phase behavior for 1/T* = 0.8, which is within the or- 
dered regime for all values of D. We see that spherical 
micelle phases are only stable in regions where tethers 
are close together (i.e. small 9) and the nanosphere di- 
ameter is "large." Conversely, lamellar phases are stable 



in regions where tethers are well spaced (i.e. large 9) 
and nanosphere diameter is "small." This diagram also 
demonstrates the strong dependence of most of the struc- 
tures on nanosphere diameter, e.g. a small decrease in 
nanosphere diameter may drive the system from ZnS/D 
to NaCl/SC or from AD to AG. The notable exception 
is the tetragonally ordered cylinder phase. Additionally, 
we can conjecture that many of these phases may be rel- 
atively insensitive to non-idealities present in experimen- 
tal systems. For example, the L/M phase forms over a 
very wide range of 9 values and thus may be relatively 
insensitive to polydispersity in bond angle. Likewise, the 
tetragonally ordered cylinder phases (TC/C and TC/T) 
form for all D values explored, thus they may be less 
sensitive to polydispersity in nanosphere diameter. 
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FIG. 8: Phase diagram for D verses 8 at 1/T* = 0.8, created 
from the data presented in figure [2] Phase boundaries are 
approximate and drawn to guide the eye. 



The progression of phases with respect to decreas- 
ing 9 is similar to changing the length of the middle 
block (B block) in a linear triblock copolymer. In linear 
triblock copolymers, as the length of the middle block 
is decreased, we typically see a change from lamellae 
[3H1 EH HE] to a tricontinuous structure (e.g. the alternat- 
ing gyroid [48] or alternating diamond phase [46, 47] [49]) 
to alternating tetragonal cylinders [46 , 48 to CsCl struc- 
tured micelles [3E1I3E]- In both cases, whether changing 
9 or the length of the middle block, the net result is that 
the A and C blocks are brought closer together, con- 
straining the possible tether configurations and resulting 
in a phase transition. Additionally, we see that decreas- 
ing the size of the nanospheres drives the system to form 
structures with less intcrfacial curvature (e.g. lamellae). 

The lamellae, alternating diamond, alternating gyroid, 
and alternating tetragonal cylinder phases are all well 
known in the linear triblock copolymer literature and 
their formation in the DTNS system is not entirely sur- 
prising; in these cases, the DTNS behave very similar to 
linear triblocks and the geometry of the nanoparticle ap- 
pears to have little impact on the overall bulk structure. 
However, linear triblock copolymers have been shown to 
form CsCl ordered spherical micelles [48] [50] rather than 
the NaCl or ZnS ordered micelles we find for DTNS. It 
has been calculated for linear triblock copolymers that 
CsCl has a lower free energy than NaCl [50]. Thus, we 
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would expect that DTNS might also form the CsCl struc- 
ture and the fact it is not present is somewhat surpris- 
ing. As we discussed in reference [T7], a nanosphere has 
no conformational degrees of freedom as compared to an 
equivalent flexible polymer with the same excluded vol- 
ume; the radius and effective volume (i.e. shape or mass 
distribution) of a nanosphere is constant, whereas the 
radius of gyration and effective volume of a flexible poly- 
mer may vary based on solvent conditions, temperature, 
or volume fraction. As a result, the flexible middle block 
in a triblock copolymer allows the A and C blocks a large 
amount of conformational freedom by changing its con- 
formation. In the DTNS system, we remove many of 
the degrees of freedom of the middle block by replacing 
it with a nanosphere that has a fixed volume/geometric 
contribution and by including a bond angle constraint, 
6. This limits the conformational entropy of the A and 
C tethers as compared to the triblock system. As such, 
these constraints lead to the stabilization of the NaCl 
structure rather than CsCl. 




FIG. 9: (a) Schematic of a star triblock copolymer, (b) 
Schematic of DTNS for 0=30°. 



Similarly, the geometric and architectural constraints 
of the DTNS are important to the formation of the 
ZnS/D structure. In the region where we find the ZnS/D 
phase, the DTNS building block closely resembles a star 
triblock copolymer, as sketched in figure [9] In reference 



[IT] we showed that the geometry of the nanosphere was 
crucial to the formation of the ZnS/D structure over the 
commonly observed [6,6,6] columnar structure. The ge- 
ometry of the nanosphere induces a larger amount of cur- 
vature than an equivalent flexible block in a star triblock 
copolymer, resulting in spherical, rather than cylindri- 
cal, morphologies |17j . This conclusion is supported by 
recent work that showed ZnS ordered spherical micelles 
for a system composed of a mixture of ABC star tri- 
block copolymers and B homopolymers |51J . The ad- 
dition of the B homopolymer acts to increase the bulk 
of the B block in the star triblock rather than just the 
length. This is similar to the impact of the geometry of 
the nanospheres and highlights the importance that the 
distribution of volume (i.e. the shape) can have on the 
overall structure. 

V. SUMMARY 

Our results show that the structural phase behav- 
ior of DTNS is a function of the directionality and 
strength of the tether-tether interactions and diameter 
of the nanosphere. We have shown that DTNS can pro- 
duce unique structural arrangements of both tethers and 
nanoparticles that are similar in bulk structure to tri- 
block copolymers. We have demonstrated a route to 
form diamond and SC networks of nanoparticles, two 
structures highly sought for photonic applications [52 
as well as sheets, cylinders, and tetragonal meshes of 
nanospheres. Overall, we have shown that the use of soft- 
matter tethers with directionality can be used to produce 
highly ordered periodic structures that would not neces- 
sarily be expected of either equivalent flexible polymer 
systems or pure nanoparticle systems in the absence of 
tethers. 
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Grant No. DE-FG02-02ER46000 and the University of 
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